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Diamagnetic Screening of the Magnetic Field in Accreting Neutron 
Stars II - The effect of polar cap widening 
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ABSTRACT 

Recently, we have proposed a model for the screening of the magnetic field of an accreting 
neutron star by the accreted material flowing from the polar regions towards the equator and 
sinking there underneath the surface (Choudhuri & Konar 2002). In this model it was assumed 
that the flow pattern remained stationary over time. However, as the surface magnetic field 
weakens, the accretion takes place over a wider region around the pole, making the flow more 
radial and isotropic. In the present work, we extend this two-dimensional model to include 
the time-dependence of the flow of the accreted material. The final radial flow is found to be 
less efficient in screening the magnetic field compared to the initial tangential flow. After an 
initial phase of rapid decay, the magnetic field slowly reaches an asymptotic value when the 
accretion becomes nearly isotropic and radial. Assuming the initial extent of the polar cap 
to be ^ S'^-IO'^, a simple geometric argument suggests that the magnetic field should decay 
by 3 — 4 orders of magnitude before stabilizing to an asymptotic value, consistent with the 
magnetic fields observed in millisecond pulsars. 
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1 INTRODUCTION 

Observations suggest a connection between the low magnetic field 
of some neutron stars (binary/millisecond pulsars) and their being 
processed in binary systems, indicating an accretion-induced field 
decay in such cases. Radio pulsar data and their statistical analy- 
ses indicate that binary as well as millisecond and globular clus- 
ter pulsars, which almost always have a binary history and shorter 
spin periods, possess much lower field strengths, down to ~ W^G 
compared to the isolated pulsars which appear to have much larger 
magnetic fields (~ lO^^G ~ W^^G) and longer spin periods (see 
Lorimer (2001) and references therein). Evidently the accreting 
material induces a rapid dissipation of the magnetic field and also 
spins the neutron star up by bringing in additional angular momen- 
tum. 

We do not have any a priori knowledge of the interior currents 
supporting the magnetic field of a neutron star. Depending on the 
generation mechanism, the field could either be supported by the 
crustal currents (Blandford, De Campli, & Konigl 1979) or by the 
Abrikosov fluxoids of the proton superconductor in the stellar core. 
The evolution of the magnetic field in these two cases are expected 
to be different. Accordingly, two classes of models have been pro- 
posed for the evolution of the magnetic field in accreting neutron 
stars (see Bhattacharya (2002), Konar & Bhattacharya (2001) for 
brief reviews) - one relating the field evolution to the spin evolu- 
tion and the other attributing the field evolution to direct effects of 
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mass accretion. Evidently, the starting point for these models is to 
assume different kinds of initial field configurations. 

Another possible method of field reduction is to screen it by 
the accreting material. As the highly conducting accreting plasma 
settles and spreads onto the surface of the neutron star, it could pro- 
duce a diamagnetic screening effect, burying the stellar field un- 
derneath it. Interestingly, this mechanism would not depend on the 
location of the field in the stellar interior and would work irrespec- 
tive of the nature of the interior currents. Though this mechanism 
was suggested quite early on (Bisnovatyi-Kogan & Romberg 1974; 
Taam & van den Heuvel 1986; Romani 1990; Romani 1995), it is 
only recently that the problem is being investigated in some quan- 
titative detail. 

It turns out that the accretion flow in the polar region of a 
magnetized neutron star is extremely sensitive to various magneto- 
hydrodynamic instabilities and it is difficult to make even a qual- 
itative assessment of the effectiveness of screening without a full 
three-dimensional computation which is yet to be attempted. A 
one-dimensional plane-parallel model by Cummings et al. (2001) 
has indicated that the diamagnetic screening is ineffective for field 
strengths above ~ 10^"G as well as for accretion rates below ~ 1% 
of the local Eddington rate and at lower accretion rates the field can 
diffuse through the accreting matter. Recently, Melatos & Phinney 
(2001) have calculated the hydromagnetic structure of a neutron 
star accreting symmetrically at both the magnetic poles as a func- 
tion of the accreted mass, starting from a polytropic sphere plus 
a centrally located magnetic dipole and have evolved the config- 
uration through a quasi-static sequence of two-dimensional Grad- 
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Shafranov equilibria with increasing accreted mass. They find that 
the accreted material spreads equator-ward under its own weight, 
compressing the magnetic field into a thin botmdary layer and 
burying it everywhere except in a narrow, equatorial belt. Ac- 
cording to their calculation the magnetic dipole moment scales 
as Bo'^M°^*Mo^ where Bo, M and Ma are the initial field 
strength, the rate of accretion and the total accreted mass. 

In an earlier work (Choudhuri & Konar (2002) - Paper I 
henceforth) we have presented a two-dimensional model to demon- 
strate the mechanism of diamagnetic screening. A neutron star with 
strong magnetic field undergoes polar-cap accretion. The material 
falling onto the polar cap flows horizontally along the surface to- 
wards the equator. The flows coming from the two polar caps meet 
near the equator and then sink below the surface. A subsurface 
poleward counter-flow is expected in the layers immediately be- 
neath the surface layer of equator-ward flow. Deeper inside, the 
material settles radially onto the core. We obtained an analytical 
expression describing this flow in Paper I and studied the evolution 
of the magnetic field subject to such a flow, kinematically, using 
a numerical code. In the absence of magnetic buoyancy, we fotmd 
that the magnetic field gets screened in the very short time scale of 
the surface flow. On the other hand, if magnetic buoyancy in the 
molten surface layers is included, the magnetic field is screened 
in the longer time scale of the interior flow. For a typical neutron 
star, these two time scales are estimated to be 1 year and 10^ years 
respectively. As magnetic buoyancy is likely to be important, we 
concluded that the magnetic field would be screened in the longer 
time scale of 10^ years. 

Observations of millisecond and binary pulsars also indicate 
that their magnetic fields undergo little or no field reduction af- 
ter the phase of active mass accretion is over. It is believed that 
the decaying magnetic field reaches an asymptotic value and then 
stops decaying any further. One of the limitations of Paper I has 
been that the calculations presented there threw no light on this as- 
pect of the problem. We used a time-independent velocity field and 
consequently the magnetic field continued to decay indefinitely. In 
reality, the velocity field is expected to be highly time dependent. 
As the magnetic field becomes weaker and the magnetic pressure 
in the polar region drops, it is possible for the accreting material 
to flow through an increasingly larger region around the pole. In 
other words, the polar cap widens with time in an accreting neu- 
tron star. Eventually, the magnetic field becomes too small to be 
able to channeUze the flow of the accreting material and the accre- 
tion becomes spherically sjmimetric. The aim of the present paper 
is to self-consistently incorporate this time dependence of the ve- 
locity field into the evolution of the magnetic field. We show that 
spherical accretion is less effective in screening a magnetic field 
compared to accretion through polar caps. Hence, as the polar cap 
widens and the accretion tends to be more isotropic, the screening 
of the magnetic field becomes less efficient and the decay of the 
magnetic field is arrested to a large extent. So, after an initial rapid 
decrease, the magnetic field reaches a phase of very slow decay 
thereafter. This result gives us a clue as to why the magnetic field 
of the neutron star attains an asymptotic value and does not decay 
any further. 

It should be mentioned here that similar conclusions were 
drawn for the case of crustal currents undergoing accretion-induced 
ohmic dissipation (Konar & Bhattacharya 1997; Konar & Bhat- 
tacharya 1999a) where a purely spherical accretion was assumed 
to be operative at all times. In that case such behaviour of the mag- 
netic field depends crucially on the nature of the detailed struc- 
tural physics of the neutron star crust. The same scenario prevails 



even for a spin-down induced expelled flux subject to ohmic dissi- 
pation in the crust (Konar & Bhattacharya 1999b). These calcula- 
tions augment the present work in which we trace the evolution of 
the magnetic field in the surface layers as the nature of the accre- 
tion changes with time. At later times, when the accretion becomes 
spherical the calculations presented in the earlier models would be- 
come relevant. 

The layout of this paper is as follows. The mathematical for- 
mulation of the problem is presented in §2. In §3, we investigate 
the evolution of the magnetic field with time-independent flow, with 
different widths of the polar cap. In §4, we obtain the time evolution 
of the surface field with a time-varying flow pattern, adjusting itself 
to the change in the field strength. However, the results obtained 
in §4 are based on an assumption of the existence of a non-radial 
flow till very deep layers of the star, which is rather unphysical and 
serves only the purpose of demonstration. Therefore, in §5 we go 
on to look at the problem where the material flow is confined to 
a thin crustal region. There has been some disagreement regard- 
ing the flow pattern, first presented in Paper I, in particular about 
the presence of a reverse flow in the deeper layers. Therefore, in 
§6 we have demonstrated that the flow velocity proposed by us is 
quite general and incorporates flow velocity without a reverse flow 
as well and have compared our results with those obtained by using 
this kind of flow pattern. Finally, we summarize our conclusions in 
§7. 



2 MATHEMATICAL FORMULATION 

We follow the mathematical formulation developed in Paper I and 

present a brief summary here for the sake of completeness. The 
magnetic field of the neutron star evolves according to the induction 
equation: 

^ = Vx(vxB)-^Vx(ivxB), (1) 
at 47r a 

where a is the electrical conductivity of the medium (see, for exam- 
ple, Choudhuri (1998), Parker (1979)). Assimring an axisymmetric 
poloidal field, allowing us to represent the magnetic field in the 
form B = V X (A(r, 9)e4,), we find that A evolves according to 
the equation: 

f + 1(v.V)(M) = ,(v^-^)a, (2) 

where 77 — /A-kg and s = rsin^. Evidently, it is the poloidal 
component of v that affects the evolution of A. We integrate equa- 
tion (2) subject to the following boundary conditions. The field 
lines from the two hemispheres should match smoothly at the equa- 
tor, requiring dA/dO = 0at6' = 7r/2. To avoid a singularity at the 
pole, we should have A = at ^ = 0. The magnetic field matches 
a potential field at the surface of the neutron star, whereas the lower 
boundary allows free advection of the field below. See Paper I for 
details of how these boundary condition are implemented. 

Assuming a polar-cap accretion, we expect the material accu- 
mulated in the polar region to cause an equator-ward flow at the 
surface in both the hemispheres. We assume this flow to be con- 
fined in a shell which is primarily in the liquid part of the crust. 
Near the equator, the flows originating from the two polar regions 
meet, turn around and sink under the surface resulting in a pole- 
ward counter-flow immediately beneath the top layer. Eventually 
this material settles radially onto the core. It should be noted that 
the pole- ward counter-flow as weU as the radially inward flow takes 
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Figure 1. Flow velocity, pv, and its divergence in a right-angular slice (0 
9 ^ 7r/2, 0.25 < r < 1.0). We have used r,„ = 0.7.5, = 0.5 for 
this picture. The panels marked la, 2a, 3a, 4a correspond to flow velocities 
withe = 0.0, 0.25, 0.5, 1.0 and flie panels marked 1&, 26, 36, 46 are the 
divergences corresponding to those respectively. 



place mainly in the solid crystalline region. The new material en- 
ters the region of interest only at the polar cap and therefore the 
flow-velocity has non-zero divergence only near the polar region. 
Apart from that V.(pv) should vanish everywhere else. 

The analytic form of a velocity field, which has all the charac- 
teristics outlined above, is given below. 
For rm <r < Vs, 




.+ 




r + 



Tb 



X (erf(V;8 cos 61) - cos^^ . 



For r < Tb, 
Ks 

' ^2 



PVr 



(5) 

) 

(6) 
(7) 

pv^ = 0. (8) 

Notice that the parameter 7 in v°{l3, 7, r^, r-m, fs) defines the size 
of the polar cap, the angular extent being given by ~ 7"^^^^. 
Whereas, the parameter /3 determines the angular width of the 
down-flow that sinks inward below the equator. The flow is pre- 
dominantly horizontal in the surface layer < < r"s. In the 
layer rb < r < r-m immediately below that, there is a pole-ward 
counter-flow. In this layer, the flow also tends to become more and 
more radial as it travels toward the deeper layers. Finally, below rb, 
the material velocity becomes completely radial, flowing inward. 

The flow in the innermost layers is nothing but the radial com- 
pression experienced by the deeper layers of the star due to an 
increase in the mass. Therefore, the magnitude of the velocity is 
related to the rate of accretion by the condition if 3 = M /An. En- 
suring continuity across r = rm and r = rj, we can relate K\ and 
if 2 to K3 and find that V.(pv°) = everywhere below r = rm- 
However, the divergence in the uppermost layer above r = rm is : 




Erf(v^cos6») ) (9) 



,(3) 



V.(pv° 



providing for a source of material only around the polar region in 
the upper layer. The panels la, 16 of Fig. 1 corresponds to this flow 
velocity and its divergence, defined above. (In Fig. 1 of Paper I 
we have shown this velocity field w°{f3 = 10.0, 7 = 10.0, = 
0.75, rb = 0.5) and its divergence V.(pv°). Please note that there 
is a typographical error in Paper I, where the values of both /3 and 
7 are mistakenly given as 1.0 instead of 10.0.). 

As the magnetic field in the polar region weakens due to 
screening, (a) the polar cap widens and (b) the velocity field tends 
to become more radial. Eventually, when the magnetic field is suffi- 
ciently weak, the accretion becomes isotropic and spherical instead 
of being polar, with the flow velocity becoming completely radial 
everywhere. We represent both of these effects through a single pa- 
rameter c taking values within the range (0, 1). A larger value of c 
implies a wider extent of the polar cap and a more radial inflow. The 
effect of the widening of the polar cap can be made self-consistent 
by relating 7 and c in the following manner: 
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ain + cAey 



(10) 



Since y^^''^ gives tlie angular width of the polar cap, we can see 
that c = corresponds to an angular width of ~ ^^min, whereas 
c = 1 corresponds to ^min + AS. In this work, we have assumed 
^min = 0.2 and A9 =1.4 such that the polar cap angle has a range 
of 10° — 90°. Through 7, the flow velocity voiP, 7, rt, rm,rs) 
too depends implicitly on c. In order to make the flow velocity more 
isotropic with the widening of the polar cap, an isotropic part is 
added to the velocity field. Therefore, the expression for the veloc- 
ity field, for a time-dependent magnetic field strength and hence a 
time-dependent polar cap area, is given by (only showing the de- 
pendence on c) 



V = (1 - c) v" (7(c)) + cv^ . 
Here v""" is the purely isotropic part given by 
K3 



pv 



-^(l-exp(5^)) 
r-' \ \ rs — rmJ J 



(11) 



(12) 



It is evident that pv""^ vanishes at the surface and reaches an asymp- 
totic value of —Ki / somewhere beneath the surface. It should be 
noted here that we are concerned with velocities inside the neutron 
star and not with those with which the accreting material may hit 
the surface of the star (which is certainly not zero at the surface). 
As freshly-accreted material makes the new surface up, the original 
surface moves deeper down. The above expressions define the ma- 
terial flow with respect to the surface at any given instant of time. 
Note that equations II and 12 imply a radial velocity of —Kz/r^ 
below r = rh, which is independent of c. Thus the rate of inflow of 
material to the neutron star core (which is determined by M and is 
related to A'a) does not change with c. 

Fig. 1 shows the velocity field and its divergence for different 
values of c. It can be seen that the polar cap (indicated by the region 
where V.(pv) is substantially different from zero, implying depo- 
sition of accreted material) becomes larger with increasing values 
of c and the flow velocity also becomes more isotropic. Finally, 
for c = 1.0, the velocity field becomes fully isotropic and radially 
inward everywhere. 

We now discuss how c at a certain instant of time depends 
on the magnetic configuration at that instant. The structure of the 
magnetosphere around an accreting neutron star is a well-studied 
subject (see, for example, Shapiro & Teukolsky (1983)). The extent 
of the polar cap area is limited by open field lines on the surface 
of the star. The last open field line is the one that goes through 
the Alfven radius, ta, in the equatorial plane. Therefore, using the 
equation for dipolar field lines, we have 

sin6lp =\j-) (13) 

where 6'p is the angle between the last open field line and the mag- 
netic axis, and therefore, the angular extent of the polar cap. At the 
Alfven radius of an accreting system, the ram pressure of the in- 
falling material equals the pressure of the magnetic field. Now, the 
ram pressure at the Alfven radius is 

Pr.u. = \p{TA)V\rA), (14) 

where p(rA) is the density and V{rA) is the velocity of the accret- 
ing material at Alfven radius. Assuming matter to be in approxi- 
mately free fall, we have V{t\\) = {2GM /rA)^^'^ ■ The density of 
accreting material at Alfven radius is p(rA) = M /'iirV{rA)r\. 
Since the field is dipolar in nature, the field strength at the Alfven 
radius is given by 



B{rA) 



(15) 



where Bs is the field strength at neutron star surface. Equating the 
ram pressure and the magnetic pressure, one obtains 



(2GM) 



(16) 



Assuming the rate of accretion (M) and other stellar parameters 
(rs , M) to be constant, we then have the following relation between 
the angular extent of the polar cap and the surface field strength: 



sin0(t)p oc Bs{ty 



-2/7 



(17) 



In our model, the velocity field automatically adjusts to the 
change in the extent of the polar cap through a variation in 7. It 
follows from equation (10) that the opening angle of the polar cap 

0{t)p, given by 7^"^^^, is 



eit)p = e,nin + C{t) AO . 



(18) 



Assuming ^min to correspond to the initial surface value of the 
magnetic field Bs{t = 0), it follows from equation (17) and equa- 
tion (18) that 



sin(0n 



c{t) AO) 



Bsjt) 
B,{t^O) 



-2/7 



(19) 



It may be noted that equation (19) is valid only during the phase 
of accretion when the opening angle increases with the weakening 
magnetic field. When c{t) equals one the accretion is completely 
spherical implying the extent of the polar cap to be 7r/2. At this 
point, the value of c.{t) is frozen and no further dependence of the 
magnetic field is included in c{t) and the flow velocity profile. 

To obtain the evolution of the magnetic field with time, we 
solve equation (2) with the velocity field described above. This ve- 
locity field, in turn, depends on the strength of the magnetic field 
through the parameter c. In most of our calculations, we have in- 
cluded magnetic buoyancy in the manner prescribed in Paper I, i.e, 
by including a radially upward velocity in the top layer. It has al- 
ready been mentioned in Paper I that the numerical program used 
for this work has been adopted from that developed in the context 
of solar MHD which formed the basis of several calculations in- 
volving solar magnetic fields (Dikpati & Choudhuri 1994; Dikpati 
& Choudhuri 1995; Choudhuri, Schussler, & Dikpati 1995; Choud- 
huri & Dikpati 1999; Nandy & Choudhuri 2001; Nandy & Choud- 
huri 2002). The details of the numerical scheme can be found in the 
Appendix of Paper I. 



3 MODIFIED FLOW VELOCITY 

Before investigating the evolution of the magnetic field for a time- 
dependent flow velocity (resulting from a time-dependent size of 
the polar cap), we look at some representative cases with steady 
flow patterns corresponding to different sizes of the polar cap. For 
this, we use the modified form of the material flow velocity given 
by equation (11). The results of this section give an indication as 
to how the variation in the flow pattern affects the evolution of the 
magnetic field. 

As in Paper I, we take the radius of the neutron star to be 
the unit of length. The unit of time is then fixed by equating the 
constant K3 in the expression of velocity to unity. Our previous es- 
timate shows that the equator ward flow at the surface takes place 
only through a narrow layer of 100 m, which is ~ 1% of the stellar 
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Figure 2. Initial field configuration, assuming a dipole confined to an outer 
region defined by 0.25 ^ r ^ 1 and 0^6^ 7r/2. 
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Figure 3. Evolution of the mid-latitude surface field with time. The curves 
a, b, c, d, e correspond to c = 0.0, 0.25, 0.5, 0.75, 1.0 respectively, ri = 
0.05 and = 50 for all the curves. 

radius (Paper I). It is difficult to handle such a narrow layer in a nu- 
merical scheme. Moreover, it has been seen that most of the qualita- 
tive behaviour can be understood quite well by using a wider flow 
layer. Consequently, we have assumed r-m = 0.75, rb = 0.5 for 
our calculations in this section and the next. Representative results 
with a narrow layer will be presented in §5. We take the diffusivity 
to be ?7 = 0.05. As shown in Paper I, the magnetic field is more 
or less frozen in the fluid for this value of 77 and the results do not 
change on decreasing r) any further. The parameter Vmh prescribing 
the magnetic buoyancy is set to the value 50.0. For all our calcula- 
tions we assume /3 to be equal to 10.0. 

We assume the initial currents to be entirely confined to the outer 
surface layers of the star and start with an initial magnetic field con- 
figuration shown in Fig. 2. This configuration is allowed to evolve 
in time with steady flow velocities corresponding to different sizes 
of the polar cap. The flow patterns for different extents of polar cap 
are obtained by taking different values of the parameter c (as de- 
picted in the velocity fields shown in Fig. 1). Fig. 3 shows how the 
magnetic field at a point on the surface at the mid-latitude decays 
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Figure 4. Field configuration at t = 1.0 starting from the initial configu- 
ration shown in Fig. 2. The panels marked a, b correspond to c = 0.0, 1.0 
respectively (corresponding to the curves a and e of Fig. 3). In both the 
cases, r} = 0.05 and v^^b = 50 are assumed. 

with time. It is clear that the decay slows down with increasing c, 
i.e. with increasing width of the polar cap. As the velocity field 
becomes more radial, it is less effective in screening the magnetic 
field. Fig. 4 shows the magnetic field configurations at time t = 1.0 
for c = and c = 1.0. It is evident that when c is zero, the effect 
of a reverse flow in the interior layers is quite prominent. Whereas, 
when c equals unity, i.e. the accretion is purely spherical, the field 
configuration shows only the effect of a radial compression. The 
fact that a more radial velocity field is less effective in screening 
the magnetic field, helps us to understand the results obtained with 
time-dependent velocity field more readily. 



4 TIME-DEPENDENT FLOW VELOCITY 

As accretion progresses, the magnetic field becomes weaker (due 
to screening, and also due to the ohmic dissipation in case of the 
crustal currents) and consequently the polar cap widens accord- 
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Figure 5. Evolution of the mid-latitude surface field with time for a time- 
dependent flow pattern (soUd curve) with r) = 0.05 and v^nb = 0. The 
dashed curve shows the evolution of c{t) with time. 



ing to equation (17). This translates to equation (19) which gives 
the time-dependence of the parameter c, in terms of the time- 
dependence of the field strength at the surface. 

Fig. 5 shows the evolution of the magnetic field at a point on 
the surface at the mid-latitude with time. It is clear that the decay 
slows down with time. As the polar cap widens out with time, the 
parameter c also increases. Since the magnetic field decay is slower 
for larger c, as we saw in the previous section, the evolution seen in 
Fig. 5 is quite expected. Fig. 6 shows the magnetic configurations 
at three instants of time, demonstrating the physical nature of the 
distortion and effective screening of the magnetic field. As expected 
from the results of the previous section, we see that gradually the 
field configuration shows the effect of a reverse interior flow less 
and less. 

One of the important characteristics of binary and millisecond 
pulsars is that their surface field strength is weaker by a factor of 
about 10^ - 10'' compared to that of the normal pulsars. Therefore, 
it has been stressed that the magnetic field strength should reduce 
by such a factor during the phase of active accretion. We propose an 
extremely simple and elegant explanation for this effect. It is seen 
that the screening becomes much weaker after the accretion has 
become isotropic. Hence, effectively, the magnetic field decreases 
until it is weak enough for the polar cap to open up completely 
to 90° and then the field reduction slows down. It follows from 
equation (19) that the ratio of the initial to the final field strength 
should be of the order of ~ (sin 90*'/ sin ^min)^''^- Then, if 6'min 
IS about 5° - 10°, we obtain a factor of of 10^ - 10'' between the 
initial and the final surface field. Therefore, if the accretion onto a 
neutron star starts with a polar cap width of about 5" — 10°, the 
decrease in the magnetic field comes out to be of the right order of 
magnitude. For the case presented in Figs. 5 and 6, we have taken 
6'min = 0.2 such that the factor sin"'^/^ 6imin is ~ 2.83 x 10^. 
Notice that in Fig. 5 the magnetic field decays approximately by 
the logarithm of this factor (~ 2.45) before the decay essentially 
stops. 

It should be noted here that the argument offered above does 
not account for the fact that the field decay continues via ohmic 
dissipation, if the currents reside in the crustal region, even after 
the accretion has become purely spherical (Konar & Bhattacharya 
1997). Detailed modeling by (Konar & Bhattacharya 1999a) has 
shown that the ohmic dissipation is effectively stalled and the sur- 
face field strength reaches an asymptotic value only when material 




Figure 6. Field configuration at intermediate times starting from the initial 
configuration shown in Fig. 2. The panels marked a, 6, c correspond to t = 
0.015, 0.05, 0.1 respectively. 
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Figure 7. Field configuration at intermediate times starting from panel (a), 
confined to a region with r ^ 0.8 for r) = 0.05 and v^\j = 50. The panels 
marked 6, c correspond to t = 0.05, 0.1 respectively. 
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Figure 8. Evolution of the mid-latitude surface field with time for a time- 
independent flow pattern (solid curve). The initial field configuration is con- 
fined to a thin region as shown in the panel (a) of Fig. 7. r) = 0.05 and 
f mb = 50 has been assumed here. The dashed curve shows the evolution 
of c(t) with time. 



of mass AM ~ lO^^M© is accreted. Therefore, if the time-scales 
of the opening of the polar cap is comparable to that of the accretion 
of AM, the above argimient would have a stronger claim to valid- 
ity. However, a more detailed calculation of the screening process, 
taking into account the micro-physics of the neutron star crust, is 
required for a definite answer in this connection. 



5 NARROW CRUSTAL REGION 

It has already been mentioned that all our previous calculations are 
done by assuming the material flow happening in a wide subsur- 
face region, with the bottom boundary at rt = 0.5, than what is 
expected in reality. This is done for the ease of niunerical calcula- 
tion. However, in order to have an idea of the field evolution in such 
a narrow region here we reproduce some of the calculations of §4 
confining the initial field configuration as well as the material flow 
in a region defined by Vb = 0.90 and rm = 0.95, starting with an 
initial configuration shown in Fig. 7. But, it should be noted that in 
a real neutron star the flow region is confined to rm =0.99 which 
is beyond the viability limit of our numerical code. 

Fig. 8 shows the evolution of the magnetic field at a point on 
the surface at the mid-latitude with time. It is clear that the decay is 
much faster compared to that in the situation where the flow hap- 
pens in a much deeper region (Fig. 3). This is due to the fact that 
the distortion of the field lines occurs in a narrower region allowing 
for faster dissipation of the field through ohmic dissipation as can 
be interpreted from Fig. 7. Moreover, the field saturates at a value 
~ 2.5 order of magnitude smaller than the original, which is about 
the same as what we see in Fig. 5. This lends credence to our ar- 
gument that the field reduction basically depends on a geometric 
factor arising out of the polar cap opening and is independent of 
the depth of flow region. By extrapolating the results of this and the 
previous section, we expect that even in the case of a very narrow 
flow region appropriate for a neutron star, the magnetic field would 
first decay rapidly by a factor similar to what we have found above 
and would then saturate to an asymptotic value. 
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Figure 9. Material velocity without a region of interior reverse flow. 
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Figure 10. Evolution of the mid-latitude surface field with time for a time- 
independent flow pattern. The curves marked a, b correspond to a material 
movement with and without (Fig. 9) a reverse flow. 77 = 0.05 and = 
has been assvmied for both the cases. 

6 VELOCITY WITHOUT COUNTERFLOW 

It needs to be mentioned here that there is no clear consensus re- 
garding the nature of the material flow. Many are of the opinion 
that there should be no counterflow and the flow in the surface layer 
should smoothly go over to a radial flow in the inner region (Gep- 
pert U., Konenkov D., Spruit H. C, private communication). The 
analytical form of the flow velocity, proposed by us, encompasses 
this possibility too. In the limit of Vm — > the flow takes this 
particular form as can be seen from Fig. 9. 

Fig. 10 shows the evolution of the magnetic field at a point 
on the surface at the mid-latitude with time. Comparing the case 
without a counterflow to the case with such a counterflow, it is ev- 
ident that the decay is slower in the first case. Since, in absence of 
a cotmterflow in the region immediately beneath the surface layer, 
the stretching and the subsequent distortion in the magnetic field 
is less, the effective decrease in the surface field strength is also 
expected to be smaller than in the case with a counterflow present. 
Moreover, as can be seen in Fig. 1 1, a flow pattern without an inter- 
mediate region of counterflow seem to give rise to a sharp bimch- 




Figure 11. Internal field configuration att = 0.1 starting from the initial 
configuration shown in Fig. 2. The panel marked a corresponds to the ve- 
locity without a reverse flow shown in Fig. 9 whereas the configuration in 
panel b has been evolved with a velocity pattern with a reverse flow (panel 
la of Fig. 1). 

ing in the internal field configuration characteristic of the region 
in which the tangential surface flow changes over to a radial flow. 
This kind of bunching may perhaps lead to a faster dissipation of 
the currents at later times which we have not investigated at present. 
To sum up, the evolution of the magnetic field is qualitatively simi- 
lar even when there is no counterflow, although the factor by which 
the magnetic field decreases is somewhat smaller. 



7 CONCLUSIONS 

Many of the questions regarding the screening of the surface mag- 
netic field of a neutron star still remain unanswered. We have tried 
to address this issue by presenting the first two-dimensional anal- 
ysis of the problem. In our previous work, we have considered a 
steady flow pattern. In the present paper, we extend this to encom- 
pass the more realistic case of the accretion flow becoming more 
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and more spherical with a reduction in the surface field strength. 
Therefore, with this work, the mathematical formulation of the 
problem can be said to be complete. However, a more detailed in- 
vestigation of this problem incorporating the internal structure and 
the micro-physics of the neutron star crust is needed, and wc expect 
to present the results of that investigation in a future communica- 
tion (Konar 2003). 

The main conclusion of this work is that the diamagnetic 
screening of the surface magnetic field becomes progressively less 
effective with time. As the strength of the surface field decreases 
due to screening, the magnetic field is less efficient in channeling 
the material flow and the accretion becomes more and more radial. 
Such radial in-fall of matter does not drag and bury the field lines 
unlike in the case of a horizontal flow. This fact is borne out by our 
calculations, and it is seen that the decay effectively stops after the 
polar cap has fully opened up. We estimate that the screening by 
accreting material should reduce the field by a geometric factor of 
about ~ 10^ — 10'' reaching an asymptotic value thereafter Hap- 
pily, these numbers match with what is observed in binary and mil- 
lisecond pulsars (compared to the field strength in normal pulsars), 
providing motivation for further investigations in this direction. 



This paper has been typeset from a TgX/ ET^ file prepared by the 
author. 
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